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ABSTRACT 


An extensive aircraft underflight program was conducted along the Skylab ground- 
path for the purpose of documenting wind, wave, and atmospheric conditions affecting 
the amplitude of the active and passive microwave signatures. The S-193‘microwave 
system senses a roughness parameter at the ocean surface that is proportional to the 
surface windspeed. The exact relationship between the wind and this roughness para- 
meter is the subject of continuing investigations. 

In the case of the active portion of the system, the intensity of off-nadir 
backscatter from the ocean is thought to be primarily determined by the amplitude 
of short gravity/capillary waves and has been shown to be strongly a function of 
azimuth relative to the surface wind direction. The passive side of the instru- 
ment senses the naturally emitted (and reflected) microwave energy and is propor- 
tional to the RMS slope and percent foam coverage of the ocean. 

NOAA, NASA, and USAF aircraft were equipped with a variety of environmental 
sensors in an attempt to specify the surface conditions affecting the satellite 
sensors as well as active and passive microwave sensors intended to calibrate the 
Skylab instrument. The aircraft program is described, and some comparisons of sat- 
ellite and aircraft results are presented. The principal result of the comparison 
of active radar is that direct inferences of the surface windspeed are possible, 
but subject to considerable scatter, and that this scatter appears to be due to in- 
teraction between long gravity and short Bragg waves and backscatter due to rain as 
well as errors in correcting for azimuth dependence. It is shown that oq for inci- 
dence angles of -50° increases both with windspeed and with increasing energy level 
of the high-frequency gravity waves that, themselves, are proportional to both the 
local wind and fetch in a manner that is not uniquely determined by the windspeed. 

An unforeseen opportunity to observe a Pacific hurricane by both Skylab and 
NOAA aircraft has contributed to the development of a simplified wave forecasting 
scheme applicable to hurricanes, and more general conditions, which combines the 
better qualities of both spectral and height/period forecasting techniques. The 
implication of this result to the SEASAT program is that quite large data inputs 
in both the time and space domain can be handled using existing computers and should 
produce a forecast of comparable or superior quality to existing spectral techniques, 
but in shorter time steps. Horizontal polarization data obtained by the aircraft in 
Hurricane Ava, and in other experiments, which led to this development are presented. 
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INTRODUCTION 


The S-193 experiment aboard Skylab was intended to test the concept of remote 
determination of surface wind conditions on a worldwide basis as an aid to improved 
environmental forecasts. A part of this experiment included an aircraft underflight 
program to calibrate and validate the inference of surface wind conditions from the 
satellite measurements of radar backscatter, ctq» ^nd upwelling naturally emitted 
microwave energy, T,. 

a 

A number of significant developments have evolved as a result of this combina- 
tion of satellite and aircraft studies of the environment, which will be touched 
upon herein. It is impossible, however, to discuss in depth all aspects included 
in the scope of this program. Results presented here are preliminary, and some con- 
clusions may be modified as additional data sets are considered. The reader is, 
therefore, asked to be tolerant of missing details and tentative conclusions that, 
hopefully, will be firmed up in future reports. 


AIRCRAFT INSTRUMENTATION 


The aircraft measurements were obtained from instrumented NOAA, NASA, and USAF 
C130 aircraft. All aircraft were equipped with a basic environmental package con- 
sisting of an inertial platform for windspeed determination, a Barnes Infrared Ra- 
diometer for sea surface temperatures, a laser wave profiler for wave measurements, 
a Cambridge Dew Point Hygrometer, a Rosemount Air Temperature Probe, and a vertical 
camera for white-cap photography. The NOAA Cl 30 was additionally equipped with a 
three- frequency passive microwave system at (14.5 GHz), X (8.35 GHz), and L (1.8 
GHz) bands. This system was used for dual polarization measurements at X and 
bands and horizontal polarization measurements at L-band. Because the NOAA aircraft 
was not available for the third launch period, a cooperative arrangement with the 
Air Force 53rd Weather Reconnaissance Squadron of the Air Weather Service was nego- 
tiated by NOAA, and the inertial platform, laser, microwave, and camera systems 
were transferred to the Air Force C130 for the SL-4 underflight program. Figure 1 
shows the NOAA Cl 30 aircraft with passive microwave antennas extended to the in- 
flight operating position. 

The NASA JSC C130 aircraft was similarly equipped except that the microwave 
system was both active and passive (time-shared) and centered only at 13.5 GHz. 

The NASA aircraft participated in the underflight program during all three launch 
periods. 


SATELLITE INSTRUMENTATION 


Skylab S-193 Instrument . The S-193 experiment aboard Skylab consists of an 
instrument capable of operation in one of three modes: (1) short pulse altimetry, 

(2) radar backscatter, and (3) a passive radiometer mode. The instrument is a single- 
frequency device centered at 13 GHz (Ku band). In the so-called RADSCAT configura- 
tion, the instrument alternately switches between the active and passive mode. It 
is capable of scanning in the along-track and across-track directions, and thus ob- 
tains dual polarization measurements of radar backscattering cross section, ap, and 
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upwelling naturally emitted microwave energy, which is directly proportional to the 
apparent, or antenna, temperature Tg, as a function of incidence angle. 

Radar Backscatter . The mechanisms responsible for determining the level of 
radar energy backscattered by the ocean surface are primarily, but not limited to, 
the amplitude of resonant Bragg waves for incidence angles between about 30° and 
80°, and the RMS slope distribution for near vertical incidence angles (Wright, 
1968; Valenzuela, 1968). The Bragg condition is defined as 


T — ” “ ^0 01 Kq cos 

water 
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where 6i and 02 are the incident and scattering angles, Kq = 2tt/X microwave. In the 
case of S-193, these waves correspond to cm wavelength ocean waves. The detailed 
behavior of Bragg waves and also their effect on backscatter has been the subject of 
considerable study in laboratory experiments in recent years (cf. Rouse and Moore, 
1972; Duncan, Keller, and Wright, 1974; Keller, Larson, and Wright, 1974; Keller and 
Wright, 1975; Reece and Shemdin, 1974; and many others). The principal result of 
these studies is that the backscattered doppler radar spectrum for a particular radar 
frequency and incidence angle is not a unique function of windspeed but rather i<^, to 
some degree, a function of fetch, duration, and presence of swell as well as wir ■ - 
speed. 


The Wave Spectrum . It has been suggested that the S-193 RADSCAT data are af- 
fected by both long and short ocean waves. In the following sections, an attempt 
will be made to evaluate S-193 data in terms of the wave spectrum. A brief descrip- 
tion of the spectrum and its behavior during growth and swell situations is there- 
fore appropriate. 


Wave conditions can be fully described by the two-dimensional energy spectrum, 
A^(f, 0 ), and integration over direction ( 0 ) and frequency, f, yields the total energy 


00 7 T 
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which is proportional to mean surface displacement, 2 a 2 . 

For a wind blowing off a shoreline, the fetch is defined simply as the distance 
from shore to the downwind measurement site. Figure 2 presents the behavior of the 
one-dimensional energy spectrum for seven fetch locations during offshore wind con- 
ditions as observed in the JONSWAP experiment (Hasselmann et. al., 1973). The in- 
set of Figure 2 depicts the meaning of the five parameters suggested by Hasselmann 
et. al. as convenient for describing the important characteristics of the spectrum. 
It can be seen from this figure that the peak frequency, tm, grows and migrates to- 
ward lower frequencies as fetch increases. Hasselman et. al. found that the 
sharpness parameter, y> and the width parameters a^, at, are essentially independent 
of fetch. The Phillips parameter, a, however, general ly decreased with increasing 
fetch. The behavior of the spectral scale parameters can also be represented 
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conveniently in nondimensional form by incorporation of the local windspeed, Uiq. and 
gravity. Thus, Figure 3 presents the behavior of nondimensional energy 

and nondimensional frequency 


versus nondimensional fetch 


Nondimensional izati on in this manner has properly accounted for the windspeed depen- 
dence of the behavior of total energy and peak frequency with fetch. The interre- 
lationship between a, E, and f^^^ is therefore apparent and was studied in detail by 
Hasselmann et. al. (1975). They proposed f^ as a convenient means of describing 
the stage of development of the wave spectrum in an average sense. Figures 4 and 

5 from this study show the behavior of E and a with respect to f and X. The lines 
denoted 
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are representative of the momentum entering the wave field (Hasselmann et. al . , 
1973), with an approximate mean value of 

C = 10"^, 

U) 

or about 20 percent of the total momentum transferred to the ocean by the wind. 
Variations about the mean are attributed to changes in the local wind conditions. 
An example of varying the wind by a factor of 1.5 for a particular f^ is included 
in the figures. 


RESULTS 


Because the principal purpose of the S-193 experiment was to remotely infer 
surface wind conditions, it is appropriate to first consider examples of the data 
from a typical pass as a prelude to more extensive analysis of the complete data 
set. Figure 6 presents the results of a pass in the Gulf of Mexico during SL-2 
and is a typical example of variable low windspeed conditions. The S-193 active 
radar backscattering cross sections, (oq), and passive microwave antenna tempera- 
ture, Tg, for horizontal polarization and an incidence angle of 50° are shown 
plotted as a function of latitude along with the surface wind conditions as de- 
termined from the NOAA aircraft tinderfl ight. 
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Considerable variability can be seen for ao> especially around 19° and 22° N. 
The data for those latitudes less than 23°, however, were not corrected for azimuth 
dependence and therefore show scatter due to changes in wind direction as well as 
speed. The higher latitudes, however, are of most interest because they were accom- 
panied by extensive ground truth. It can be seen that the surface winds decreased 
significantly between 24° and 26° N, while oq decreased only slightly between 24° 
and 25.5° N and actually increased at 26° N. The horizontally polarized antenna 
temperature agrees well with aircraft determinations and also increases signifi- 
cantly at 26° N. Because of the presence of many rain showers observed by the air- 
craft to be in the area of 26° N, these increases are tentatively attributed to 
this factor. This pass, therefore, implies significant corrections would normally 
be required to account for the presence of rain. 

During SL-2, a unique opportunity to observe high wind conditions developed 
with the appearance of a hurricane in the eastern Pacific southwest of Acapulco, 
Mexico (Ross et. al., 1974), As plans were being formulated for observing the 
storm with the S-193 system, the NOAA C130 aircraft deployed on the 6th of June to 
Acapulco, refueled, and flew a 7.5-hour mission into the storm. Figure 7 is a NOAA- 
2 composite satellite view of Hurricane Ava showing the flight track of Skylab as it 
conducted a data pass with S-193 operating in the side-looking solar inertial mode. 
Unfortunately, the storm was rapidly moving away from the subsatellite track and it 
was not possible to obtain measurements in the region of maximum winds. Figure 8 
is an example of active and passive measurements of Ava obtained at incidence angles 
of 42.5° and 50.5° along with estimates of surface winds obtained during the air- 
craft penetration and rainfall rates estimated from the NIMBUS-5 satellite 19.35 
passive microwave system (Wilheit, 1972). It can be seen from this figure that 
there is an increase in oq and T^ in the higher wind areas of the storm with the 
highest antenna temperature occurring in the zone of heavy rainfall. 

High wind conditions were also obtained during the third launch period. Figure 
9 shows the subsatellite track for data obtained on the 9th of January at incidence 
angles of 0° and ± =50° along with the NOAA surface analysis for 1800Z. This situa- 
tion is particularly interesting because the windspeed varied from 7.5 to 30 m/sec 
from the beginning to the end of the sampling period with little rainfall reported 
except in the region of the front. Figure 10 presents the variation in oq and Tg 
for +47.6° and -50.5° incidence angle along with surface winds estimated from an 
isotach analysis based mainly on ship reports and the NASA JSC aircraft measure- 
ments. A qualitative comparison of these data sets, together with that from Hurri- 
cane Ava and the June 11 pass in the Gulf of Mexico, strongly suggests a first- 
order dependency on surface wind conditions but with scatter. 

Unfortunately, due to damage to the S-193 antenna occurring during an SL-3 
extravehicular excursion, the antenna pattern was altered in an unknown fashion. 

It is therefore risky to include the 9 January data set with those obtained during 
SL-2 when calculating windspeed dependency. This case was therefore treated sepa- 
rately and is shown plotted against windspeed in Figure 11. oq and Tg data ob- 
tained during SL-2 are shown in Figure 12 plotted against windspeed. These data 
sets are restricted to those data passes described above, which were accompanied 
by an aircraft underflight. Aircraft-determined antenna temperatures included 
show good agreement with S-193. The surface winds attached to each satellite data 
point are judged to be accurate to about 1 to 3 m/sec in the case of SL-2 and 3 to 
5 m/sec in the case of the SL-4 pass of 9 January. The judgment of accuracy in- 
cludes the effects of mesoscale variability in the local wind conditions and inac- 
curacies associated with aircraft-determined winds and ship reports in the case of 
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the 9 January data. Unfortunately, no study is known to the author that compares 
ship reports (mostly visual estimates) to continuously recorded and calibrated 
winds averaged over 10 to 30 minutes. Therefore, one is left with a judgment. 


DISCUSSION 


The above data sets strongly suggest a useful first-order relationship between 
surface winds for both ao and Tg subject to some degree of scatter. There are sev- 
eral sources of errors in both parameters contributing to the scatter that crops up 
in the processing. In addition, there may be errors due to invalid atmospheric as- 
sumptions (or errors in the corrections required for a particular assumption), ran- 
dom errors in inferring the 10-meter windspeed from aircraft measurements, as well 
as natural variability in local wind conditions. The above errors inherent to this 
data set probably cannot be reduced any further. 

One potential source of error that can be addressed, however, lies in the lack 
of uniqueness of the active or passive signature due to the variability in possible 
wave conditions that may be present for a given windspeed. Reece and Shemdin (1974), 
in a study conducted in a wave tank, showed that the high-frequency waves, for a par- 
ticular fetch, are windspeed dependent, but the absolute energy level for a particu- 
lar windspeed is reduced with the addition of a low-frequency component (swell) and 
that the amount of the reduction is proportional to the amplitude of the low- 
frequency component. Mitsuyasu (1971) showed similar reduction in the high- 
frequency gravity region of the wave spectrum with the introduction of swell. 
Hasselmann et. al . (1973), showed that the Phillips constant (a) (Phillips, 1958), 
which determines the energy level of the f"^ region of the wave spectrum, decreases 
with increasing nondimensional fetch as discussed earlier. In order to assess the 
possible importance of the gravity wave spectrum in this data set, it is desirable 
to consider ag as a function of some observable parameter of the wave field that 
varies with fetch in a well-behaved and predictable manner. 

The nondimensional peak frequency, was shown earlier to be a particularly 
useful parameter that well describes the stage of development of the wave,spectrum. 
Aircraft measurements of the wave spectrum were used directly to specify f^ for 
Skylab data sets obtained on 5 and 11 June 1973. The data for Hurricane Ava, how- 
ever, present a special problem because the aircraft measurements were not obtained 
at the exact subsatellite point. In order to specify the peak frequency, it was 
desirable to develop some technique of estimating these parameters from the air- 
craft data set. The hurricane wind fields are circular in nature, however, and 
the fetch relationship needed to infer is ambiguous and arbitrary. Furthermore, 
the position of the aircraft measurements relative to the eye of the hurricane were 
concentrated in the rear quadrant, whereas most of the satellite positions were to 
the right of the hurricane center. Figure 13, however, presents some of the wave 
data obtained, plotted along with spectra from the North Sea measured under similar 
wind conditions, but fetch limited, and in the North Atlantic during fully developed 
conditions. The resemblance between the spectra is striking. The most significant 
feature of the Ava spectra, however, is the general lack of "swell," which would ap- 
pear as a secondary peak in the spectrum. This is especially significant when the 
unidirectional assumption required in processing the aircraft data is considered. 
Because we assume all waves are moving in a direction parallel to the aircraft 
flight track, swell from some other direction is moved toward lower frequencies in 
the mapping process to fixed coordinates, leading to an unrealistic "broad" 
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appearance to the spectrum with multiple peaks. Hurricane Ava was a superhurricane 
by any criteria and had a record low pressure for eastern Pacific storms of 914 mil- 
libars. Fortunately, another such storm of which wave measurements have recently 
become available was Hurricane Camille, one of the worst hurricanes to ever strike 
the coast of the United States as it went ashore near Mobile, Alabama, in August of 
1969. Camille's eye dimensions, maximum winds, forward velocity, and central pres- 
sure were virtually identical to those of Hurricane Ava. Measurements of wave con- 
ditions in relatively deep water were obtained by a consortium of oil companies. 

Some of the data were reported recently (Patterson, 1974, and Hamilton and Ward, 
1974) and were used in this study. Figure 14 presents the relative positions of 
Camille and Ava wave data to the eye and wind field of Ava as determined from the 
NOAA aircraft flight. It is fortunate that Camille was so similar to Ava as the 
windspeeds measured on the oil company platform were biased low because of poor 
anemometer exposure for a storm approaching in the direction of Camille. Ava's 
winds were therefore used to specify the Camille wind field, which avoids the dif- 
ficulty of introducing another empirically based technique to specify 10-meter ane- 
mometer winds. The wave data from Camille and Ava were then nondimensional ized and 
plotted against the nondimensional radial distance from the eye and are presented 
in Figure 15. It can be seen that a simple power law reasonably well describes the 

radial behavior of both E and f^. Nondimensional peak frequencies for the Ava data 

set were therefore calculated from the expression. 


where 


and r is the particular radial distance from the subsatellite point to the eye of 
Ava. The Ava data set was then combined with data from 5 and 11 JuneJ973 and is 
shown in Figure 16. It can be seen that oq varies considerably with f„,. It can be 
argued that, because the windspeed is included in the calculation of f„,, it is dif- 
ficult to correctly separate the windspeed dependency from stage of development. 

To aid in this separation of dependencies, a multiple regression analysis was per- 
formed according to the equation 


Z = ao + a^x + &2y 

letting Z = oq, x = Uiq, and y = f^. The constants ag, ai,_and a .2 were found to be 
-29.5, 0.20, and 59.4, respectively. Such a dependency on fm is much greater than 
expected or would be predicted on the basis of wave tank experiments cited earlier. 
The same multiple regression analysis was performed on the 9 January data set, based 
on hindcast fm, for 47.6° incidence angle, and yielded values of -15.7, 0.39, and 
17.5 for the same constants. The hindcast performed assumed limited fetch, but un- 
limited duration, and the rather low wave heights reported suggest an underestimate 
for the values of fm for the higher winds. Also possible is a bias in the case of 
the hurricane data set due to backscatter from rain. The 9 January data set, how- 
ever, seems to confirm a fetch dependency and both sets considered together suggest 
that controlled high-wind, variable-fetch experiments should be performed to accur- 
ately infer high winds from measurements of oq. 
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A similar treatment of the fetch effect in the case of the passive microwave 
measurements is also indicated and is underway. 


RELATION TO SEASAT 


The success of the simplified approach suggested by Hasselmann et. al. (1975) 
for specifying the evolution of the wave spectrum depends on successful parameteri- 
zation of the nonlinear interactions that control the exchange of energy within the 
spectrum that, in turn, are very sensitive to local gustiness in the surface wind 
conditions. However, because it is not necessary to deal with the entire directional 
spectrum for each grid point in a numerical forecast scheme, it will be possible to 
increase the density of grid points and decrease the time steps involved in fore- 
casting waves. For example, a typical spectral model consisting of 17 frequencies 
and 15 directions for 512 grid points in the North Atlantic requires 130 000 storage 
locations, whereas the simplified approach, expanded to account for two swell sys- 
tems, can increase the grid density to 5012 and require only 30 072 locations. 

Such a forecasting approach is needed for rapid assimilation of satellite data and 

will likely be in operation by the time the SEASAT satellite is launched. At this 

time, the First Global GARP Experiment will also be underway and provide a unique 
opportunity to test the SEASAT concept. 


CONCLUSIONS 


The Skylab S-193 experiment has proved that active and passive microwave sen- 
sors can be used to infer surface winds but are subject to scatter and a decreasing 
sensitivity with increasing windspeed in the case of the active radar, and bias due 
to rainfall with little sensitivity to lower windspeeds in the case of passive micro- 
wave signatures. These results, therefore, suggest that combined active and passive 
systems with a weighted averaging process (employing polarization dependencies) being 
used to infer the local wind might reduce some of the scatter due to random errors 
and should be tested with Skylab data. The results further suggest that a parameter- 
ization of the wave spectrum may be necessary in order to further reduce the scatter 
in oq. The Skylab data set contains most of the data needed to test these hypotheses 
and could lead to satellite determination of both the windspeed and surface wave 
spectrum by judicious use of active and multi frequency passive microwave systems. 
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Figure 1. The NOAA C130 aircraft show passive microwave sensors in the extended position 
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Figure 2. Growth of wave spectra for offshore wind conditions during JONSWAP-69. Fetch 
increases from Station 5 through Station 11. Three shape parameters 0 |^, 

and y, and the scale parameters f^^ and a suggested by Hasselmann are shown in 

the inset, y is simply the ratio of the energy at f^ to that which would be 

predicted by the Pierson-Moskowitz (1964) form of the spectrum. 



Figure 3. Dimensionless parameters E and versus dimensionless fetch X for high wind- 

speeds (15 to 25 m/sec) as determined by aircraft experiments in the North Sea 
and off Cape Fear, N. C., on 27 January (Ross and Cardone, 1974, and Barnett 
and Wilkerson, 1969). 
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Figure 4. The behavior of E and during growth. A sudden change of windspeed by a 

factor of 1.5 results in a departure from the mean followed by readjustment to 
mean condition along the curved lines indicated. Nondimensional fetch relation 

ships X = E = u^- 2 " particular f^ are included. The upper and lower bound- 

“3 -5 

aries (C^ = 10 and = 10 ) represent the limits of momentum being trans- 

ferred to the wave spectrum. 





Figure 5. The Phillips parameter a versus Corresponding fetches 5 are indicated by 


the tick marks. The behavior of a during readjustment to a change in local 
wind (U) conditions is along the curved lines. 
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Figure 7. Hurricane Ava, June 6, 1973, NOAA-2 composite showing track of Skylab satellite. 
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SKYLAB SURFACE CHART JANUARY 9, 1974 

Figure 9. Surface meteorological analysis for 9 January 1974 pass. Subsatellite positions 
of S-193 under study are shown for the indicated incidence angles. 
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Figure 10. S-193 data during the 9 January pass along with surface winds determined from 

the surface analysis. 
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Figure 11. Radar backscatter, oq versus windspeed for 9 January. 



(qp) (X.) 



U,o (m/d) 


Figure 12. Radar backscatter and antenna temperature versus windspeed during SL-2. 
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Figure 13. Wave spectra measured during the aircraft penetration of Hurricane Ava along 
with North Sea spectra for fetch-limited 25 m/sec winds. 
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Figure 14. Surface wind analysis for hurricane Ava constructed from aircraft data along 

with relative positions of wave data obtained in Ava and for Hurricane Camille. 



E, 

tut 


and peak frequency, f^, versus radial 
s the region of maximum winds. 







1936 



Figure 16. Radar backscatter, oq, versi^s stage of wave development as represented by non- 
dimensional peak frequency for the SL-2 data set. 



